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ABSTRACT 

The results of an experimental calibration of the NASA Glenn 
Research Center 16” Mass-Flow Plug (MFP) are presented and 
compared to a previously obtained calibration of a 15” Mass-Flow 
Plug. An ASME low-beta, long-radius nozzle was used as the 
calibration reference. The discharge coefficient for the ASME nozzle 
was obtained by numerically simulating the flow through the nozzle 
from the WIND-US code. The results showed agreement between the 
15” and 16” MFPs for area ratios (MFP to pipe area ratio) greater than 
0.6 but deviate at area ratios below this value for reasons that are not 
fully understood. A general uncertainty analysis was also performed 
and indicates that large uncertainties in the calibration are present for 
low MFP area ratios. 

INTRODUCTION 

Large-scale model supersonic aircraft inlets are routinely tested in 
NASA Glenn Research Center’s (GRC) larger wind tunnels such as the 
8x6ft and 10x10ft Supersonic Wind Tunnels (SWT). Two key 
objectives of such tests are to ensure that the air supplied by an inlet 
matches the mass-flow demand of an engine and that the inlet operates 
in a stable condition over a wide range of mass-flows. Thus, accurate 
measurement of inlet mass-flow is of primary importance. Due to its 
relatively simple mechanical design, the mass flow through an inlet is 
often controlled and measured by a calibrated conical mass-flow plug 
(MFP). For larger models, NASA GRC uses a 16” diameter MFP 
which was fabricated by modifying a pre-existing 15” MFP. 
Calibration of the 16” MFP plug was performed in-house using an 
ASME low-beta, long-radius nozzle as the reference measurement. 
The procedure and resulting calibration, as well as a comparison to the 
calibration of the previous 15” MFP are reported herein. 

NOMENCLATURE 

A = Area (in 2 ) 

AR = MFP area to cold-pipe area ratio (Eq. 20) 

C D = Discharge coefficient 


cl = ASME nozzle diameter (in) 

D = Diameter (in) 

g c = Proportionality constant (32.17405 lb- ft/(lbf- s 2 ) 

Lj = ASME nozzle throat length (in. Figure 3) 

M = Mach number 

MFP = Mass Flow Plug 

p = Pressure (psi) 

R = Radius (in) 

Rj = ASME nozzle ellipse major radius (in. Figure 3) 

R 2 = ASME nozzle ellipse minor radius (in. Figure 3) 

Rail - Gas constant for air (53.3525 ftTbf/(lb-R) 

Re rf = ASME nozzle Reynolds number 

R e ds = Scaled ASME nozzle Reynolds number 

Re d s = Re rf X 10E-06 

s = Slant height of the frustum of a cone (in, Figure 5) 

T = Temperature (R) or throat tap location (in. Figure 3) 

U = Velocity (ft/sec) 

w = mass flow rate (lb/s) 

x,y,z = Cartesian coordinates 

x,r, 6 = Cylindrical coordinates 

x P ] ug = Measured position of MFP (in) 

x 'piug = Effective position of MFP (in. Figure 5 and Eq. 23) 


Greek Symbols 

P = Ratio of ASME nozzle-to-approach pipe diameter 

y = Ratio of specific heats 

8Xj = Uncertainty of measurand X, 

p = Density (lb/ft 3 ) 

<j) = MFP cone half-angle (deg, Figure 5) 


Subscripts 

a 

amb 

8 * 

i 

n 


= Pertaining to the average component 
= Pertaining to ambient conditions 
= Pertaining to boundary-layer blockage (displacement 
thickness) 

= Pertaining to ideal conditions 
= Pertaining to the normal component 
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noz = Pertaining to the ASME nozzle 

pipe = Pertaining to the cold-pipe 

plug = Pertaining to the MFP 

pot = Pertaining to potential flow conditions 

t = Pertaining to total conditions 

v = Pertaining to the vertical component 

Superscripts 

* = Pertaining to choked flow conditions 


Pamb(^) 

w- 4 * 



CALIBRATION PROCEDURE 

The calibration of both the caged 15” and 16” mass-flow plugs 
were carried out with the experimental setup shown in Figure 1. The 
four primary components of the calibration rig are the American 
Society of Mechanical Engineers (ASME) inlet flow nozzle assembly, 
an adapter spool piece, the cold-pipe assembly, and the mass-flow plug 
assembly. The entire calibration assembly was mounted vertically on 
the floor of NASA Glenn’s lOxlO-ft Supersonic Wind Tunnel and 
protruded through a sealed opening in the ceiling. The ASME nozzle 
was fed by atmospheric ambient pressure and temperature air while the 
wind tunnel exhausters were used as a vacuum source to establish a 
flow through the assembly. A photo of the ASME nozzle on the top of 


the 10 x 10-ft. SWT is shown in Figure 2 and a schematic of the 
geometry is shown in Figure 3. A photo of the 16” MFP attached to the 
cold-pipe is shown in Figure 4 and a schematic of the plug is shown in 
Figure 5. 

The spool piece is an adapter to match the ASME nozzle exit 
diameter to the cold-pipe diameter under test. The cold-pipe is simply 
a 96” long constant-area pipe whose primary purpose is to allow for a 
flow-settling length between the inlet exit and the mass-flow plug. 
Near the forward end is a spherical-nosed centerbody supported by 
three equally spaced struts. The centerbody nose is intended to 
simulate the fan spinner on a turbojet engine. Attached to the forward 
end of the cold-pipe is a spool piece that houses eight compressor face 
rakes. The rakes lie in the plane of the tip of the centerbody nose. The 
inside edge at the exit plane of the cold-pipe is sharp to fix the position 
of the mass-flow plug choke point. 



Figure 2. ASME nozzle assembly. 



A picture of the 16-inch mass-flow plug attached to the cold-pipe 
is shown in Figure 4. The cone is actuated by a hydraulic cylinder and 
concentricity of the cone with the cold-pipe is maintained by wedge 
bearings on three of the six knife-edged support struts (cage). A 
position transducer (potentiometer) can be seen on the top of the 
assembly in Figure 4. 
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Figure 4. Caged 16” Mass-Flow Plug. 



Figure 5. Schematic for MFP sonic (choked) area. 


The ultimate goal of the calibration is to determine the discharge 
coefficient for the 16” MFP over an operating range typical of the 
supersonic inlet testing. The discharge coefficient for the mass-flow 
plug is defined as: 


Eq. 1 


c 


D.plug 


W 

_ noz 


w 


plug.i 


where W no _ is the actual mass flow through the cold-pipe and W* ptugi 

is the ideal inviscid choked flow through the mass-flow plug. For the 
purpose of calibrating the mass-flow plug, the actual mass flow was 
measured with an ASME Long-Radius Flow Nozzle. A description of 
this nozzle is as follows. 


ASME Nozzle 


Geometry: An ASME Long-Radius Flow Nozzle with throat 
taps was used to provide a reference mass-flow measure for calibrating 
both the 15” and 16” mass-flow plugs. The geometry of the nozzle is 
shown in Figure 3. This design conforms to the Low-P Nozzle with 
Throat Taps illustrated in Fig. II-III-14 of Ref. 2 with the following 
two exceptions. First, no approach pipe exists before the nozzle 
(D=co), hence J3=d/D=0. And second, the nozzle exit flow does not 
exhaust into a sudden expansion but rather into a shallow angle conical 


diffuser. The nozzle throat-to-conical diffuser and conical diffuser-to- 
constant area end section are both transitioned by 27.96” radii. 1 

Citing the desire to maintain fully turbulent flow in the nozzle, the 
calibration tests were performed both with and without a boundary- 
layer trip in the inlet region of the nozzle. Details of the trip were not 
supplied by Iek [1] but inspection of the ASME nozzle indicated the 
trip location and width are as shown in Figure 3. Most of the tests were 
performed with an 80 grit trip, but limited testing was also performed 
with a 36 grit trip with essentially the same results. 

Flow Equations: The flow through the ASME nozzle is 
calculated by the following equation: 


Eq. 2 


where 
Eq. 3 


W — ( ' I V’ 

noz D,noz noz,i 


w = p U -A 

noz,i / noz noz n 


and where p no . and t/ are evaluated at the throat of the ASME 


nozzle. Substitution of the ideal gas law, Mach number definition, and 
isentropic and adiabatic relations for pressure and temperature, 
respectively: 


Eq. 4 


P noz 


144 -P noz 

D . T 

v air J noz 


Eq. 5 U = M Jv • g ■ R -T 

“ noz noz v / <-> c air noz 


Eq. 6 


P noz 


P t,noz 


Eq. 7 


0 + }/ 2 ] - M Iz) 


into the expression for VV Q and rearranging yields: 
Eq. 8 


/■8c Pi, noz ' M noz ' A., 


Kir 


FT 

V 1 t,r 


1 + 


-<r+ 1) 

Hr-D 




Substituting the values for y,g c and R ajr and replacing p t 
and T t no _ with p amh and T amb , respectively, gives the final form: 


Eq. 9 


=0.91884- 


P amb 


■ M noz- A «oz 




+ 0.2 -M: 


Note that this diffuser geometry differs from the description given by 
Iek [1] where the diffuser is described as being a “cubic contour.” 
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where M is the Mach number of the flow at the throat of the 

noz 

nozzle and is given by the isentropic relation: 


Eq. 10 


M noz = 


5- 

( 

f \ 

P amb 

f— 1 ^ 

(3.5 J 

-1 


V Pnoz J 

l 


1 

K 


y 


and p noz is the static pressure measured at the throat of the nozzle. All 
of the variables in Eq. 9 and Eq. 10 are directly measured or calculated 
from directly measured quantities. The only remaining parameter 
required to calculate the nozzle flow rate is the nozzle discharge 
coefficient (Cd moz ). 


ASME Nozzle Discharge Coefficient 


Historical Correlation: The discharge coefficient for a 
nozzle with throat taps is a function of the potential flow non- 
uniformity and the boundary-layer displacement thickness resulting 
from viscous effects. As discussed by Rivas and Shapiro [3], if the 
length of the cylindrical sections before and after the throat taps are on 
the order of or greater than 'A throat diameter, then the effects of 
potential flow non-uniformity are negligibly small. 

In the absence of a calibration, ASME recommends reading the 
discharge coefficient for a low-(3 long-radius nozzle with throat taps 
from Fig. II-III-19 of Ref. [2], The source of this curve is not 
referenced but it appears to be the “Calculated Curve” in Fig. 22 of 
Cotton and Westcott [4], No uncertainty is given for using this curve 
but based on similar analysis, it is likely to be on the order of ±1.5- 
2.0%. Much of the data used to generate this curve was based on 
incompressible flow data. In 1973, Smith and Matz [5] investigated an 
8-inch iow-P, long-radius nozzle under incompressible and 
compressible conditions and reported a Mach number dependence. 

Two-Code Computational Approach: For the March 1996 
calibration of the 15” caged mass-flow plug [1], due to the uncertain 
effect of the /? parameter and the nozzle exit diffuser on the discharge 
coefficient, the ASME correlation was not used. Instead, the 
theoretical method first performed by Smith and Matz [5] and later by 
Lahti and Hamed [6] was used to establish a discharge coefficient. 
With this method, the discharge coefficient is assumed to depend on 
two components, an inviscid part ( Co.poi ) that accounts for the radial 
pressure variation in the throat, and a viscous part (C D l j») that accounts 
for the effective boundary-layer blockage. Computer programs were 
used to calculate each component and the discharge coefficient of the 
nozzle was determined from the product of the components: 


Eq. 11 


c =c ■ c 

y -‘D,noz,1C D.pot '-D,S* 


This approach will hereafter be referred to as the “Two-Code” 
(2C) nozzle calibration. In Lahti and Hamed’s [6] study, the potential 
solution was obtained with the Stream Tube Curvature (STC) 
computer code [7], chosen because of their many years of experience 
with it and because it is driven by global mass conservation as part of 
its solution procedure making it a natural choice for mass-flow 
calibration studies. The boundary-layer correction was obtained using 
Harris and Blanchard’s code [8,9] . This code was chosen based on its 
wide use and thorough documentation. Lahti and Hamed [6] report 
that comparisons with experimentally determined discharge 
coefficients indicate that an accuracy of ±0.25% can be achieved with 
their method. 


For the calibration of the NASA Glenn ASME nozzle, Iek [1] 
followed a similar approach but used different computer codes. For 
the potential solution, Stockman and Farrell’s code [10] was used and 
for the viscous solution. Herring’s PL2 code [11] was used. With these 
codes, the theoretical discharge coefficient for the ASME nozzle was 
determined at six discrete throat Reynolds numbers. 


Table 1. Discharge Coefficient vs. Reynolds Number 


Re^ v =Re^x 10‘ 6 

Cd.koz 

1.6828 

0.98553 

3.0386 

0.98667 

3.7659 

0.98728 

4.0302 

0.98728 

4.3062 

0.98755 

4.6056 

0.98796 


For data reduction during the mass-flow plug calibrations, the data 
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in Table 1 were interpolated with a third-order polynomial function: 

c C D ^ 2C =0.982982 + 0.002081573Re rf>s - 

Eq. 12 

0.0004063900 Re^ + 0.00004062582 Re^ 

and extrapolated by setting points outside the calibration range to the 
end points of the calibration. 

RANS Computational Approach: More recently, a different 
computational approach was taken for obtaining the discharge 
coefficient of the ASME nozzle as well as the 16” MFR Rather than 
use a potential code coupled with a boundary-layer calculation, 
WIND-US, which is a general purpose 3-D Reynolds-averaged 
Navier-Stokes (RANS) code, was used for the numerical studies [12]. 
Versions 3 Alpha and 2 of WIND-US calculated the CFD derived 
discharge coefficients. WIND-US is a mature, multi-zone, structured- 
grid, compressible flow solver offering a variety of turbulence models. 
All simulations were run steady-state using the Menter Shear Stress 
Transport (SST) turbulence model [13]. The SST turbulence model 
was chosen for all of the simulations because previous studies [14,15] 
have shown that this particular model performs well for jet flows. 

An axisymmetric structured grid, composed of roughly 119,000 
grid points and divided into 17 zones, was used to model the test rig 
assembly shown in Figure 1. Grids were created for 19 different plug 
positions, ranging from x plug =2A” to 16.4” with X plug defined as 

shown in Figure 5. Zones upstream of the experimental boundary-layer 
trip location (see Figure 3) assumed laminar flow while zones 
downstream of the trip location assumed turbulent flow. An example 
grid is shown in Figure 6 while Figure 7 illustrates the solid computer 
model. 

The simulations used a constant Courant-Friedrichs-Lewy (CFL) 
number of 0.5 rather than a constant time step. To reduce 
computational costs and convergence time, parallel processing along 
with three levels of grid sequencing were used to obtain converged 
solutions. First, a coarse grid sequence included only every fourth 
point in each direction. Second, a medium grid sequence was used that 
utilized every other point in each direction. Third, a fine sequence used 
every grid point. Solution convergence was achieved when the 
computational equation residuals were reduced by several orders of 
magnitude and the mass flow rates in the bellmouth entrance and near 


Ref. 1 indicates a fourth-order polynomial was used but the actual data 
reduction program uses the third-order fit shown. 
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the cold pipe exit each varied by less than 0.09% per iteration. 
Convergence was typically obtained in a little over three hours. 



Figure 6. Sample primary grid layout, laminar zones are in 
blue, turbulent zones are in red. 



k 

TOPWORK Camera TOP 

Figure 7. 90° sector of the axisymmetric solid model. 


Initially, the WIND-US code was run with the as-tested ambient 
inlet conditions (p om *=l4.3 psia), but then to better resolve Reynolds 
number effects, additional cases were run at inlet pressures of 5.0 and 
50.0 psia. For all pressures, the ambient total temperature was set to 
520.0°R with an initial inflow Mach number set to 0.001. Whereas for 
the 5.0 psia case, the back-pressure at the MFP exit was set at 1.0 psia, 
for the two higher pressure cases, the back-pressure was set at the as- 
tested pressure of 2.0 psia. 

After the CFD analyses were completed, the data was post- 
processed to acquire the discharge coefficients for the bellmouth 
nozzle. First, the mass flow rate through the cold pipe was calculated 
from the average mass flow rate between the downstream tangential 
edge of the radial nose of the conical plug to the cold pipe exit: 

N 

w . 

pipe,x,i 


Eq. 13 


W 


pipe 


i=l 


N 


where 


Ec l' 14 w ptpvc =jjp-U x -r-dr-d0 

0 r 

This was compared to the average mass flow rate at the throat of 
the bellmouth to check for continuity. The solution was converged 
when the average mass flow rates in the bellmouth throat and just 
upstream of the cold pipe exit varied by less than 0.09% of each other. 
The discharge coefficient could then be calculated from: 


Eq. 15 C 

M D,noz 

where W noTj is given by Eq. 9. 



Red (X10 6 ) 

Figure 8. ASME nozzle discharge coefficient. 



0.1 1.0 10.0 

Re d (X10 6 ) 

Figure 9. ASME nozzle discharge coefficient, constant Mach 

number. 


The data from the WIND-US computations are plotted in Figure 8. 
For reference, the ASME curve, the data from Smith and Matz, and the 
numerical data from Iek are also shown. For the WIND-US data, only 
points where the throat Mach number was between 0.2 and 0.8 are 
plotted. Due to the use of a boundary-layer trip on the nozzle, the 
distribution does not display the transitional behavior of Smith and 
Matz [5] data and indeed seems to agree well with the turbulent 
portion of their “Critical” curve (Re rf „ oz > 2.0E+06). Another 
observation of the current data is the saw-tooth nature of the combined 
distributions indicating a small Mach number dependence. To illustrate 
this dependence more clearly, data at constant Mach numbers of 0.3, 
0.5 and 0.7 were extracted and are plotted in Figure 9. Although the 
dependence on Mach number is small for the current data, the trend is 
opposite that of Smith and Matz [5], but since their low Mach number 
data was in the transition region, nothing firm can be drawn from this. 
Comparing the present WIND-US data to the previous distribution 
established for the same nozzle configuration by Iek [1] shows that the 
present data trends about 0.2% higher than Iek’s. 

In order to apply easily the present data to the 16” MFP 
calibration, a correlation was developed that fits the data and accounts 
for the small Mach number dependence. This was accomplished by 
first curve fitting the M=0.5 data which resulted in the following 
expression: 
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f\ ft 

Eq. 16 c =0 993003 - 

vRe rf 

An empirical correction to this nominal discharge coefficient was 
established such that the Mach number dependant discharge 
coefficient is defined by: 


as the ratio of the actual mass-flow to the ideal theoretical mass-flow 
(see Eq. (1)). The ideal mass-flow is calculated from the following 
equation: 

Eq. 19 w* , ■ — p -U • AR ' ■ A i 

n plug,i r' pipe pipe A plug 


^9- 17 C D ,noz,WU —C d,M =0.5 ' f (M) 

where 

Eq. 18 /(M) = 0.999514 •(l-0.2M,^)' 0 0086341 ' 


where 
Eq. 20 


AR = 


A 

^ plug 

A 

pipe 


Since it is unlikely that the current ASME nozzle would ever be 
operated at anything but ambient pressure conditions, the coefficients 
in Eq. 18 were determined from the 14.3 psia data only. A comparison 
of Eq. 17 and Eq. 16 with the numerical data is shown in Figure 10a 
and Figure 10b, respectively. 




Figure 10. Comparison of ASME discharge coefficient 
correlation with numerical data. 


and where A p i ug is the choked-flow area between the cold-pipe exit 
plane edge and the conical plug. With reference to Figure 5, the 
definition of the choked area (A plllg ) of the MFP is based on the area of 
the frustum of a right circular cone where the slant height ( s a ) is 
defined as the average of the vertical distance from the cold-pipe edge 
to the MFP and the slant height of a cone whose slant height is normal 
to the MFP surface. 

Eq. 21 A, = n • S -(r + 7? ) 

“ plug a \ pipe a) 


where (see Figure 5) 


K = x' p , ug ■ Sin(^) 

Eq. 22 Sy = x'plug ' tan(^) 

S a =(s n +s v )/ 2.0 

R c, = R pipe - ■ cos (<^ - arccosO„ / S a )) 

and 

Eq- 23 x' plug = x plug + (R pipe - R plug ) / tan(^) 

where x' p i ug is the distance from the end of the cold-pipe to the 
theoretical intersection line of the plug conical surface and the cylinder 
of the cold-pipe internal diameter (see Figure 5). The choice of s a to 
approximate the sonic line (curved line in Figure 5) is arbitrary but 
must be used consistently for the calibration to be valid. 

The area ratio can be expressed in terms of the Mach number in 
the pipe by the following isentropic relation: 

Eq. 24 


Substituting Eq. 24 and Eq. 4 through Eq. 7 (with the subscript 
“noz” replaced with “pipe”) into Eq. 19 yields: 


A pi“g ( 7 + 1 ^ 

A pipe V 2 J 


7+1 

2(H) 


■M. 


1 + ^--M 2 

pipe 


-(/+!) 

2 ( 7 - 1 ) 




16” Mass-Flow Plug 

Geometry: With reference to Figure 5, the plug itself is a 
^=16.5 degree half-angle cone with a 2.18 inch nose radius. The 
downstream end of the cone is /? p; „ ; =15.800 inches in diameter 
followed by a 0.610 inch wide cylindrical shoulder. 

Flow Equations: The discharge coefficient for the mass-flow 
plug is assumed to be a function only of the plug position and defined 


Eq. 25 


W. 


i,pl“g 


rsc 

I R..:, 


f \ 

Pt,pipe 


^ „ 


■A 


7 + 1 n 


-(7+1) 

2 ( 7 - 1 ) 


plug 


The total pressure in the cold-pipe (Pt, p i pe ) is calculated from Eq. 6 
(with the subscript “noz” replaced with “pipe”). The static pressure in 
the cold-pipe (p pipe ) is directly measured with static pressure taps and 
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the Mach number in the cold-pipe ( M pipe ) is calculated from Eq. 24 
which assumes that the flow at the mass-flow plug exit is choked. 

Instrumentation: With reference to Figure 1, the test setup 
was instrumented as follows. The ambient air pressure and temperature 
were each recorded with four transducers which were then averaged to 
obtain p amb and T amb . There were a total of six throat taps on the ASME 
nozzle and eight static taps at the mid-point of the cold pipe. These 
were also averaged to obtain p mz and p pipe . The temperatures were 
measured with Type-E thermocouples and the pressures were 
measured with ESP 8400R 0-15psia modules. The data was recorded 
on the NASA GRC ESCORT data acquisition system. 

Calibration Runs: For the caged 16” mass-flow plug 
calibration runs, there were six sweeps of the mass flow plug. The first 
three were performed with the boundary-layer trip installed on the 
ASME nozzle and the last three were performed with the trip removed. 

Uncertainty Considerations 

A general uncertainty analysis was performed for both the 15” and 
16” MFP calibrations. The uncertainty in the MFP discharge 
coefficients was estimated by: 


Eq ' 26 


8X, 


where X t are the measurands which are summarized with their 
assumed uncertainty in Table 2. 


Table 2 Measurand Uncertainty 


i 

Xi 

SXi (15716” MFP) 

Units 

i 

Xplug 

0.015 

in 

2 

T 

1 amb 

2.5 

“R 

3 

Pamb 

0.0060/0.0045 

psi 

4 

Pnoz 

0.0049/0.0037 

psi 

5 

Ppipe 

0.0042/0.0032 

psi 

6 

d 

0.0025 

in 

7 

Dpipe 

0.0075 

in 

8 

Dplug 

0.0050 

in 

9 

<t> 

0.05 

deg 

10 

Cd.hoz 

0.001 

N/A 


boundary-layer trip and fully turbulent downstream. In the experiment, 
there is a finite width of grit material. These effects are not easily 
quantifiable, thus the uncertainty specified for the nozzle discharge 
coefficient is to be considered a ‘‘best guess” by the authors. Finally, 
although the ambient temperature is measured and used in the 
calculations, it does not contribute to the MFP discharge coefficient 
uncertainty. 

RESULTS AND DISCUSSION 
Caged 15” MFP Calibration 

The 15” MFP was calibrated over a range from x p/ug =0.25” to 

17.0” with and without distortion screens installed in the cold-pipe. A 
correlation was established based on the data without distortion 
screens. The discharge coefficient was assumed to be only a function 
of the choked flow area ratio (Eq. 20) and is expressed as a segmented 
third-order polynomial function where the coefficient values are given 
in Table 3: 

Eq. 27 C o,piug,2C = c o c i AR c 2 ’ c 3 ' 


Table 3. Coefficient values for caged 15" mass-flow plug. 



AR < 0.203943 

0.203943 <AR < 0.811515 

0.811515 <AR 

Co 

1.0340732 

0.96163863 

1.1845804 

Cl 

-0.58958996 

0.0047836495 

-0.60737883 

C 2 

1.1527497 

-0.030427921 

0.42665471 

C 3 

0 

0.050854378 

0 


This original calibration (Eq. 27) was based on the two-code 
ASME nozzle calibration. To update the calibration to use the more 
recent WIND-US based nozzle calibration, two approaches are 
possible. First, the data calibration data points could be recalculated 
with the WIND-US based nozzle calculation and a new fit established, 
or second, an adjustment could be applied to Eq. 27 by multiplying by 
the ratio of nozzle discharge coefficients: 

Eq. 28 C D,plug,WU — C D ,plug,2C ' ( Cd.iwz.WU I C D.IIOZ.2C ) 

The nozzle discharge coefficients based on both methods are 
plotted in Figure 11 as a function of the 15” MFP area ratio. Also 
plotted is the ratio of the coefficients which can be approximated by 
the following empirical relation: 


The measurands P amb , P noz and P pipe are actually calculated averages of 
four, six and eight individual pressure measurements, respectively, and 
the uncertainties in Table 2 reflect the measurement redundancy. In 
addition, different pressure uncertainty levels were used for the 15” 
and 16” MFP calibrations to reflect an upgraded pressure measurement 
system for the 16” MFP. The measurement uncertainty for the 
individual modules was 0.012 and 0.009psi for the 15” and 16” MFP, 
respectively [16], The uncertainty in the nozzle discharge coefficient is 
comprised of two components. The first is related to the grid and 
solution convergence and is generally small (<0.003%). The second is 
related to how well the CFD models the flow physics. In particular, 
since the nozzle discharge coefficient is primarily a function of the 
boundary-layer growth, the turbulence model will have an influence on 
the result. Further, the difference between how the boundary layer 
transitions on the nozzle will also affect the discharge coefficient. In 
the CFD, the flow is assumed laminar until the leading edge of the 


=1.0014446 + 0.017657109 • AR - 
r 

D,noz,2C 

Eq ' 29 0.00069852959/AR - 0.046818586 • AR 2 + 

2.3567944E - 05/AR 2 + 0.049185331 • AR 3 - 
4.4263225E - 07/AR 3 - 0.018204318 • AR 4 

The original and adjusted 15” MFP calibrations along with the 
data they were derived from are shown in Figure 12. The nozzle Mach 
number and Reynolds number are also plotted for reference. With 
regard to the original calibration, there are a couple of things to note 
about this distribution. At area ratios below approximately 0.2, the 
curve has a large negative slope and at the lowest area ratio, the 
discharge coefficient exceeds unity with large data scatter. At the 
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largest area ratio, the slope and scatter also increases. During the test, 
the ratio of the cold-pipe base pressure to pipe total pressure is 
monitored to ensure the plug is operating under choked conditions. 
The two highest area ratio positions are approaching flow conditions 
that challenge the assumption of a well-defined choked flow. 
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Figure 11. ASME Nozzle discharge coefficient as a function 
of 15” MFP area ratio. 



Figure 12. Original and adjusted 15” MFP discharge 
coefficient. 


The adjusted 15” MFP discharge coefficient curve fit is re-plotted 
in Figure 1 3 with the uncertainty bounds from the general uncertainty 
analysis. Also plotted is the percent uncertainty based on the assumed 
measurand uncertainties (Table 2). 
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Figure 13. 15” MFP discharge coefficient uncertainty. 


These results indicate that the measurement uncertainty exceeds 
0.5% for MFP area ratios below 0.5. Based on this calibration, the inlet 
should be sized to operate within the area ratio range 0.2<A/?<0.8 
(1.85”<x p ,„ s <10.4”). 

The relative contribution from the individual measurands is plotted 
in Figure 14. Also plotted for reference is the absolute uncertainty. The 
large uncertainty at low area ratios is primarily due to the low Mach 
number in the nozzle (see Figure 12) which is calculated from Eq. 10. 
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Figure 14. Measurand contribution to overall 15” MFP 
discharge coefficient uncertainty. 


Caged 16” MFP Calibration 

The results of the calibration of the 16” MFP for both tripped and 
un-tripped nozzle boundary layer are shown Figure 15. The nozzle 
Mach number and Reynolds number are also plotted for reference. At a 
given area ratio, due to its larger size, the 16” MFP has a larger flow 
rate than the 15” MFP. At large area ratios, this results in choking of 
the nozzle for the 16” MFP. Thus data beyond an area ratio of 0.7 will 
be excluded from further analysis as the MFP would be operating 
unchoked. 
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Figure 15. 16” MFP discharge coefficient data. 

Comparing the 16” MFP data with and without the boundary-layer 
trip shows agreement above an area ratio of about 0.5, but deviates 
significantly below this value. At an area ratio of 0.5, the nozzle 
Reynolds number is about 3.0E+06. From Figure 9, a Reynolds 
number of 3.0E+06 is consistent with the end-of-transition implied by 
the data of Smith and Matz [5] and the ASME curve [2]. Smith and 
Matz data show a large increase in discharge coefficient below this 
transition value. Since a fully turbulent nozzle discharge coefficient 
was used in the 16” MFP calibration, which was artificially low for the 
case with sweeps with no boundary-layer trip, the resulting MFP 




calibration is under predicted. In fact, what the transitional nozzle 
discharge coefficient would have been could be estimated by 
correcting the no-trip data to agree with the tripped data. 

A comparison between the adjusted 15” MFP and the tripped 16” 
MFP calibration data is shown in Figure 16. The 16” MFP calibration 
data above an area ratio of 0.7 has been removed (ASME nozzle 
choked) and the remainder of the data was fit with a second order 
polynomial: 

Eq. 30 

C Dplug = 0.943044 + 0.026915 • AR + 0.017376 • AR 2 

Also included in this plot are the uncertainty bounds from the general 
uncertainty analysis. The experimentally determined discharge 
coefficients for the 15” and 16” MFP tend to agree well above an area 
ratio of 0.6, but deviate below this ratio. In particular, the 16” MFP 
curve does not exhibit the abrupt increase in flow coefficient below an 
area ratio of 0.2. Inasmuch as the calibration curves lay at the edges of 
the others uncertainty bounds, this suggests that the differences in 
calibration are not fully attributable to measurement uncertainties. For 
both MFPs, the method of calibration produces large uncertainties at 
low area ratios. This is primarily due to the low Mach number in the 
ASME nozzle. The use of a delta pressure transducer at low Mach 
numbers would mitigate the large uncertainties significantly. 
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Figure 16. Comparison of 15” and 16” MFP calibrations. 

As previously mentioned, the numerical study with WIND-US 
included a numerically based determination of the 16” MFP discharge 
coefficient (see Figure 6 and Figure 7). This study is still in progress 
and will be reported in a future publication, but the preliminary results 
indicate a trend very similar to the experimental 16” MFP results but 
shifted up about 0 . 25 %. However, results for area ratios below 0.2 
have not yet been obtained so no comment can be made regarding the 
tail-up behavior of the 15” MFP. 

CONCLUDING REMARKS 

The results of a calibration procedure for the NASA GRC 16” 
MFP have been presented and compared to a previous calibration of a 
15” MFP. The 16” MFP was fabricated from the 15” MFP so all future 
testing will be with the 16” MFP. For the 16” MFP, a calibration has 
been established, but a difference in character between the 15” MFP 
and the 16” MFP for area ratios below 0.6 is observed. Unfortunately, 
this region of disagreement coincides with large uncertainty levels. 
The recommended path forward is to first complete the numerical 
calibration of the 16” MFP; in particular obtaining calibration data 
below an area ratio of 0.2. And second, the experimental calibration of 


the 16” MFP should be repeated with a focus on minimizing 
measurement uncertainties. 
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